Ixodes scapularis density and Borrelia burgdorferi prevalence along a residential-woodland gradient in a region of emerging Lyme disease risk

The environmental risk of Lyme disease, defined by the density of Ixodes scapularis ticks and their prevalence of Borrelia burgdorferi infection, is increasing across the Ottawa, Ontario region, making this a unique location to explore the factors associated with environmental risk along a residential-woodland gradient. In this study, we collected I. scapularis ticks and trapped Peromyscus spp. mice, tested both for tick-borne pathogens, and monitored the intensity of foraging activity by deer in residential, woodland, and residential-woodland interface zones of four neighbourhoods. We constructed mixed-effect models to test for site-specific characteristics associated with densities of questing nymphal and adult ticks and the infection prevalence of nymphal and adult ticks. Compared to residential zones, we found a strong increasing gradient in tick density from interface to woodland zones, with 4 and 15 times as many nymphal ticks, respectively. Infection prevalence of nymphs and adults together was 15 to 24 times greater in non-residential zone habitats. Ecological site characteristics, including soil moisture, leaf litter depth, and understory density, were associated with variations in nymphal density and their infection prevalence. Our results suggest that high environmental risk bordering residential areas poses a concern for human-tick encounters, highlighting the need for targeted disease prevention.

detect on a human host, posing the most significant threat to public health 48 .Moreover, the detection of multiple life stages over consecutive years provides evidence of established tick populations, which pose greater environmental risk 49 .Our findings can help municipal and public health officials inform residents in areas of increased peridomestic risk of necessary Lyme disease prevention and intervention measures.

Descriptive statistics
We completed 2736 quadrat observations from twelve grouped neighbourhood zones (three zones in each of N1-N4) in western Ottawa (Fig. 1) through nineteen rounds of tick drag sampling in 2020 and 2021.A total of 537 I. scapularis ticks were collected, of which 123 were nymphs, 300 were adults, and 114 were larvae (Table 1).More than half (54%) of the total ticks were collected in N2, three-quarters of which were from the woodland zone.The residential zone in N1 was the only area where we found no ticks of any stage over the entire study period.Also, no nymphal ticks were collected from the residential zone of N4.At least one adult tick was collected in the residential zone of all neighbourhoods except N1.Nymphs were present in the interface and woodland zones of all four neighbourhoods, with the total number ranging from 3 (N4) to 10 (N1) in interface zones and from 7 (N4) to 61 (N2) in woodlands.The difference in the number of adult ticks collected from the interface (22)  and woodland (28) zones of N1 was small, while adult ticks collected from the other neighbourhoods' interface zones were less than half of their woodland totals.Nymphs were collected most often between May and July, while adults were active between May and June and from September to the end of the study season in both years (see Supplementary Fig. S1 online).We collected larvae in all but one (N4) of the three study neighbourhoods.Though we observed none at residential sites, larval counts in N2 and N3 were higher in the interface (5 larvae collected in each) and woodlands (21 and 22, respectively).The interface zone of N1 proved to be the exception.In this zone, we found more than half (61) of all larvae collected-50 of these in one sampling effort.
Of the 123 nymphs collected, 27 (21.9%)tested positive for B. burgdorferi, while 107 (35.7%) of the 300 adults collected were infected (Table 1).We collected infected nymphal ticks from all zones except for the residential of N1 and N4, and the interface of N1.The infection prevalence among all tested ticks from interface zones was equal to or greater than that of woodland zones.Across all four neighbourhoods and both sampling years, the nymphal infection prevalence was 22% while the combined adult-nymph infection prevalence was 32% (Table 1).
Peromyscus species mice were present in all neighbourhoods and zones.In total, we sampled 812 mice of the target species.Peromyscus leucopus (i.e., white-footed mice) comprised 69% of sampled mice across all neighbourhoods, with P. maniculatus (i.e., deer mice) representing 27%.Peromyscus spp.mice were collected in the greatest abundance in N1 and N4, with 40% and 35% of the total, respectively (Table 2).We collected most mice from the residential zone of each neighbourhood, except for N3 where we collected the most from the woodland zone (Table 2).The prevalence of B. burgdorferi infection among all sampled mice was 11% (Table 2), and 96% of infected mice were white-footed mice; only 4 sampled deer mice tested positive (Supplementary Table S1).Neighbourhood prevalence of B. burgdorferi infection among all mice ranged from 6% (N4) to 18% (N2).The residential zone of N1 was the only residential area with mice that tested positive for B. burgdorferi.Among P. leucopus mice, B. burgdorferi prevalence was between 17 and 32% in interface zones and 19% and 53% in woodland zones.Overall, 57% of the infected mice we sampled were collected in woodland zones (Supplementary Table S2).The mean deer foraging intensity was 0.5 events (i.e., photos of unique deer) per 100 trap-hours, with a range from 0.2 (N4) to 1.1 (N1).All zones of N1 exhibited the highest deer foraging intensity during the period our trail cameras were in place.This neighbourhood's residential zone saw more unique events (1.8 per 100-trap hours) than any other study zone (Table 2).
We detected significant associations between all ecological characteristics and both sampling neighbourhoods and zones (p < 0.001).The most common tree species in all neighbourhoods were maples, though the difference in proportions of sampling quadrats dominated by maples compared to coniferous species was smaller in N3 and N4 (Table 3).Quadrats in interface and woodland zones were predominantly occupied by maple trees, while the dominant species among residential quadrats was more varied.Full canopy cover (69%) and understory density (41%) most characterized quadrats in woodland zones, though this was also true of N2 in general (Table 3).Shallow litter depth was most common in interface and woodland quadrats.Dry soil humidity characterized the majority of the neighbourhood and zonal quadrats.

Statistical models
Compared to residential zones, nymphal density was between 4 (interface, 95% confidence interval (95% CI): 1.19-15.03,p < 0.001) and 15 (woodland, 95% CI: 4.40-49.85,p < 0.001) times greater when controlling for other factors (Table 4), presenting a strong gradient effect across zonal habitat types.The association between nymphal infection prevalence and woodland zones was similarly strong (odds ratio (OR): 13.62, 95% CI: 2.07-89.p = 0.007), though the relationship with interface zones was not statistically significant (Table 4).Analyses that assessed the density and infection prevalence of all potentially infectious ticks (i.e., nymphs and adults) again demonstrated a strong and significant zonal gradient effect (Table 4).The density of nymph and adult ticks was nearly 22 times greater in interface zones and over 31 times greater in woodland zones relative to those in residential zones.Moreover, the infection prevalence of all tested ticks was nearly 15 times greater in interface zones than in the residential (OR: 14.90, 95% CI: 3.15-61.15,p < 0.001) and just over 24 times greater in woodlands (OR: 24.23, 95% CI: 6.37-125.14, p < 0.001).Though a substantial overlap between the confidence intervals for interface and woodlands zones existed for all metrics, their estimates relative to residential zones were statistically significant (p < 0.001).
Tick densities predicted by the average marginal effects demonstrate zonal and neighbourhood associations with the environmental risk of Lyme disease when site-specific and ecological characteristics are held at their mean or referent value (Fig. 2; Supplementary Table S3).Heterogeneous tick densities exist between the studied neighbourhoods; zones in N1 and N2 exhibited greater densities compared to those in N3 and N4 (Fig. 2A,C).Overall, predicted nymphal density was 0.05 per 100 m 2 in interface zones compared to 0.17 per 100 m 2 in woodlands.This zonal gradient on nymphal density was most evident in N2 (0.07 per 100 m 2 and 0.24 per 100 m 2 in the interface and woodland zones, respectively).The same gradient effect across zonal habitat types was present in predictions of the density of all potentially infectious ticks (i.e., nymphs and adults), both within and independent of neighbourhoods.Diagnostic plots of residuals indicated good model fits given the chosen response distributions, covariates adequately explaining the outcomes, and the absence of influential outliers (Supplementary Figs.S2-S5).

Discussion
Our study provides a fine-scale examination of Lyme disease environmental risk along a landscape gradient, demonstrating elevated risk in interface and woodland zones relative to residential paths and trails.The density and infection prevalence of nymphal ticks in these ecotonal zones were higher than those found in residential zones, highlighting the significant role they hold in the Lyme disease risk present in neighbourhoods featuring residential areas with integrated greenspace.Studies of habitat types within urban woodland parks of England found the highest densities of I. ricinus nymphal ticks in woodland edge sites, with evidence of established populations in grassland, woodland, and woodland edge habitats of urban and peri-urban greenspaces 28,29 .Conversely, in New York state, Horobik et al. found that densities of I. scapularis ticks and infected ticks were both higher in forests than at forest-field edge, concluding that edges do not pose a higher risk than the forests themselves 50 .However, each of these studies examined the entomological risk present at the edge of forests with open grassland.We examined the forest edge that specifically interfaces with the boundary of residential properties.One hypothesized mechanism underlying the effect that human-driven landscape change has on Lyme disease risk is termed "ecological release", wherein populations of a species increase due to changes that remove constraints limiting their expansion.A relevant example of this is the creation of new forest edges or clearings-specifically when these border human residential areas that provide new resources and nesting opportunities favoured by, for example, white-footed mice 22 .By accounting for possible differences in the abundance of white-footed mice in the sampled zones, which is believed to be the primary reservoir in this region 11,12 , as well as characteristics of the microhabitat known to impact blacklegged tick survival, we shed light on the role of interface zones on Lyme disease risk.Direct comparison of total tick density between the unique zones and neighbourhoods indicates heterogeneity at this scale.Densities of nymphs and all potentially infectious ticks in N2 reflected previous estimates for the Ottawa municipal region 51 .The densities identified in each of the three remaining study neighbourhoods were lower.Still, despite a small sample size of infected nymphal ticks from these three neighbourhoods, we detected an infection prevalence that also supports estimates reported previously for this region, indicating the presence of Lyme disease environmental risk 51 .It is notable, however, that the prevalence of infection among the collected nymphs was consistently higher in interface zones compared to woodlands, apart from N1 where no infected nymphs were sampled.We also note that, while evaluating the density of nymph and adult ticks together is common, doing so carries potential issues in the assessment of infection rates given that prevalence in adult ticks is typically greater, having potentially fed twice on reservoir hosts.Still, the zonal differences we detected underscore the importance of considering the ecological relationship with reservoir hosts along the environmental gradient we examined.1.
Our results from sampled Peromyscus-species mice indicated higher abundance in residential zones while B. burgdorferi prevalence was highest in woodland habitats.This supports previous evidence demonstrating greater white-footed mouse abundance at sites with less forested area within a 500-m radius yet lower B. burgdorferi prevalence than their within-woodland kin 44 .Together, these results are demonstrative of the differences in scale-based associations of landscape edge effects.Though the difference in nymphal infection prevalence we measured for the interface zone was not statistically significant, the findings of Piedmonte et al. corroborate a similarity in prevalence among all potentially infectious ticks through their exploration of multi-scale habitat effects on tick-borne pathogen prevalence 45 .Continued surveillance in residential-woodland ecotonal habitats, where the potential for human-tick encounters is high, will improve our understanding of this relationship.And, while the B. burgdorferi prevalence of white-footed mice observed in our data appears to support its role in enzootic transmission within this region, the involvement of other potential host species remains understudied.Additional rigorous evaluation of host community composition would help to more fully characterize the transmission dynamics between these zonal settings 52,53 .
Several limitations of this study should be considered along with the results described above.Weak relationships between some ecological observations and the evaluated outcomes may be due to the generalization of each 50 m 2 sampling area to discrete observations.We cannot dismiss the possibility that environmental measurements were subject to observer bias that resulted in misclassification.While it is a strength of this study design to investigate microhabitat characteristics of the sampling area where ticks were collected, the precision of effects for ecological characteristics could be further improved by quantifying the same measures.For example, obtaining imagery from drones flown along sampling transects and using machine learning image processing methods to inventory and calculate distributions of tree species would be a more efficient and less costly approach to quantify these areas of risk and guide field sampling.In terms of zonal factors, the final models we evaluated were adjusted for the involvement of white-footed mice and deer in the outcome measures, but we did not account for the structure or diversity of the broader vertebrate host community; we were not able to focus on other small or medium-sized mammal species.While we adjusted for study neighbourhood in our models, differences in landscape configuration of zones within neighbourhoods that were not measured may impact the results and should be an area for further research.We also determined tick density and infection prevalence using drag sampling which is less resource-intensive than the collection of ticks from small mammal hosts but may be less efficient [54][55][56] , though its effectiveness as a method to identify I. scapularis populations in this region has been shown 57 .We followed rigorous standardized protocols for drag sampling, but stage-differential or non-detection of ticks (false negatives) are a recognized limitation of this method.I. scapularis adults are generally easier to collect by drag sampling as some types of terrain more easily dislodge nymphal ticks (N.H. Ogden, personal communication, February 2, 2024).Finally, the prevalence calculations for some of the neighbourhood zones were based on small totals of collected nymphs and adults, which can result in wider confidence intervals due to uncertainty in the estimates 58 .Future research should continue active surveillance in these and additional neighbourhoods and other Canadian regions with emerging tick populations to further validate these estimates and their relationships.
Urban expansion involves the encroachment of residential developments into natural areas.As suburbs and peri-urban regions continue to change through this practice, it is crucial to understand how new ecotones contribute to local Lyme disease risk.Several past studies connect the quantity of wooded space in and around the residential yard to the increased risk of encounters with host-seeking ticks 35,46 .Yet these findings focus on Lyme disease risk where established residential properties border wooded areas and do not assess where new landscape change integrates existing forest with new residences, trails and pathways.Ours is among the first to quantify the effect of land cover types in different stages of urban development on local tickborne disease risk, specifically using a "One Health" approach across the residential-to-woodland ecotonal gradient.In this study, we focused on the fine-scale localization of Lyme disease risk within neighbourhoods of a region experiencing the rapid emergence of blacklegged tick populations.By simultaneously modelling the local impact of wildlife host availability, microhabitat characteristics, and stage of land development on the density and infection status of I. scapularis ticks, we demonstrated that the environmental hazard of tick-borne diseases is significant at the intersection of urban greenspace and residential properties.These results may help to inform urban planning guidelines regarding the selection and development of sites through the characterization of environmental Lyme disease risk prior to the development of new residential neighbourhoods.Continued surveillance, informing targeted control measures 59 , is necessary to ensure that local Lyme disease risk is minimized as new ecotonal "interface" areas allow human tick encounter opportunities.Finally, these effects should be investigated in other settings experiencing blacklegged tick emergence, to better understand the environmental characteristics of Lyme disease hazard and incorporate intervention and control measures into expansion plans.

Study location
Four neighbourhoods across the rural and suburban areas of western Ottawa (Fig. 1) were selected to assess tick populations, characterize the environmental factors and estimate the abundance of species involved in the B. burgdorferi transmission cycle (i.e., Peromyscus-species mice and white-tailed deer).Guided by the results of recent active tick surveillance that established the presence of black-legged tick populations 51 , we engaged in partnership and consultation with the National Capital Commission, Ottawa Public Health, and the City of Ottawa's Planning Division to arrive at an informed selection of these sites.To identify characteristics related to varying degrees of proximity to urban and residential development, we subdivided each neighbourhood into three sampling zones: residential, woodland, and interface (Fig. 1).Residential zones consisted of pathways and trails within neighbourhoods between the back yards of residential properties, while woodland zones were defined as forested areas adjacent to the developed areas of neighbourhoods and connected to residential areas by trail networks.Interface zones consisted of trails and wooded edges immediately bordering residential properties, such that homes or fences remained visible from just inside the tree line.

Tick drag sampling
Our study area included 36 sites evenly divided across the four neighbourhoods (N1-N4), with three sampling zones in each neighbourhood (i.e., residential, woodland, and interface), and three sites in each sampling zone.
Field collectors experienced in collecting ticks by drag sampling visited each site monthly from May through October 2020 and biweekly across the same period in 2021.At each site, team members followed a standard tickdragging sampling protocol, wherein a 1-m white flannel sheet is dragged along the ground to collect questing ticks 60 .The starting location for each tick drag transect was selected based on observations from preliminary site visits that a continuous linear drag remained in the same zone type without interruption and allowed for measurement from three consecutive sites across both sampling years.Field collectors stopped every 25 m to check and remove ticks from the flannel and record the geographical coordinates of the sampling location using a Garmin eTrex 20 GPS, covering 200 m.In total, the area dragged per sampling visit was 600 m 2 in each zone and 1800 m 2 in each neighbourhood.All collected ticks were returned to the lab in specimen tubes, where species and stage were identified by microscope according to standard taxonomic keys [61][62][63] .

Ecological characteristics
To ascertain any distinguishing characteristics that might influence the suitability of distinct locations for tick survival, we examined the fine-scale aspects of the microhabitat at each surveillance site.Field collectors recorded site-specific ecological measurements (Table 5) derived from several studies of blacklegged tick abundance at every second sampling location (50 m) in June of both study years.Mid-season values for all variables were used to capture and adjust for general ecological differences across all sites in our analyses.
We identified all categorical variables according to previously established classification keys 9,[64][65][66] .Field researchers also measured the percentage of canopy cover overhead at sampling locations using a spherical densiometer, according to established field protocols 67 .The same study personnel measured litter layer depth with a small ruler perpendicular to the ground, rounding to the nearest centimetre, and measured soil moisture in percentage by inserting a moisture meter probe approximately 4-6 cm into the soil.In analyses, we also treated litter depth as a categorical variable with the thresholds used by Talbot et al. 9 : none, shallow (< 2 cm), or moderate to deep (3-10 cm).

Small mammal sampling
Our trapping effort was based on a target sample size of up to 50 individuals per species (i.e., Peromyscus leucopus, or white-footed mice, and Peromyscus maniculatus, or deer mice) per zone each year, assuming 1 mouse for approximately every 10 trap-nights with additional trap-nights to account for trap failures.Trapping sessions took place in two rounds of two consecutive nights in each neighbourhood zone during the peak period of mouse activity (June to August).Ethical approval for small-mammal research was obtained from the University of Ottawa Animal Care Committee (permit ME-3079).All animal trapping and handling methods were performed in accordance with the Canadian Council on Animal Care's guidelines involving wildlife 68 .
During the two-night sampling period in each neighbourhood zone, field collectors laid 150 Sherman traps, baited with an apple cube, a pinch of oats, and a cotton round, on both sides of the trail.Team members ensured trap placement in two parallel lines 10 m apart, with 5-10 m of distance between each trap in the same line.The exact distance between traps varied depending on the habitat and availability of forested space in each location.Traps were set each evening and checked the following morning.Therefore, our sampling effort was 600 trapnights per neighbourhood zone in each sampling year, totalling 7200 trap-nights across all neighbourhoods.Deciduous trees provide heavier leaf litter and have been associated with higher levels of risk.Coniferous trees have been shown to discourage tick establishment due to the absence of leaf litter in some regions 9,74 Canopy cover % A dense or closed canopy cover can protect from weather extremes (cold, wet) and provide more leaf litter 11,15 Understory density Bare (none present; reference) Sparse (visible but no significant coverage) Dense (thick, difficult to traverse) Full (healthy, covering most of the forest floor) Habitats with a dense understory may provide better refuges for ticks from extreme weather and a grassy substrate on which to quest 9,66,[75][76][77] Litter layer depth None (0 cm; reference) Shallow (< 2 cm) Moderate to deep (3 to 10 cm) Leaf litter provides a refuge from cold in winter and keeps the microhabitat humid, which contributes to tick survival by protecting them from desiccation 9,11,15,66,77 Soil moisture % In microhabitats, soil moisture is directly related to the relative humidity, which contributes to the ability of ticks to survive 9,66,74,77 Vol:.( 1234567890 www.nature.com/scientificreports/Small mammals captured during the night were briefly removed from traps into a cotton bag to safely measure their weight before removing them for inspection.For all trapped animals of the target species, sex, reproductive condition, length and weight were recorded, attached ticks were collected, and one ear punch biopsy was taken from each ear 69 before the animal was released.Individuals for which the species could not be identified based on field guides were referred to as "P.sp." in field records 70 .For any non-target animals trapped (e.g., chipmunks, voles), we identified the species, collected any attached ticks and released them.We did not sample these diurnal species as the capture methods we employed are tailored to nocturnal animals.

Pathogen testing
We tested all biopsy ear punches, adult, and nymphal ticks for Borrelia burgdorferi s.s.using quantitative polymerase chain reaction (qPCR) assays according to previously published protocols 71,72 .In short, we extracted total genomic DNA using the QIAamp mini kit (QIAGEN, Mississauga, ON, Canada).With this, we identified Borrelia species using a duplex qPCR assay targeting the 23S ribosomal RNA.We further confirmed the presence of Borrelia burgdorferi s.s. by targeting the ospA gene.We performed amplifications with the BioRad CFX96 Real-Time PCR Detection System and sample analysis with CFX Maestro Software version 2.3 (https:// www.bio-rad.com/ en-ca/ produ ct/ cfx-maest ro-softw are-for-cfx-real-time-pcr-instr uments, Bio-Rad Laboratories, Hercules, CA, USA).

Deer density
As a proxy of white-tailed deer density in our models of tick outcomes, we calculated deer foraging intensity in each zone as captured by trail cameras.We installed one trail camera at each neighbourhood zone (12 total locations) at the end of June each year to collect images of large mammals active in the area.Every two weeks until the middle of October, field collectors relocated each camera 200 m along the tick drag sampling site, resulting in 8 camera trap sites per neighbourhood zone.We filtered all collected photographs to select only images of white-tailed deer and totalled the number of use events such that a unique deer captured by the camera counted as one event.To avoid overestimation due to double-counting, photographs with deer captured within a 10-min interval of the previous capture were considered a single event to allow for the individual's activity to cross in and out of the camera's detection zone.We estimated the foraging intensity at each camera location by dividing the number of events by the total number of camera-hours at the site, as per previously documented calculation methods 73 .We then averaged this number across all twelve installation locations for the camera to arrive at a "foraging intensity" estimate for deer activity in each neighbourhood zone.

Statistical analysis
We linked ecological data to the total number of I. scapularis ticks of each stage from the current and preceding sampling locations.This resulted in 50-m-long rectangular quadrats as the unit of analysis, with four quadrats in each of the three sampling sites of all neighbourhood zones (144 in total).Each quadrat was dragged for questing ticks 19 times across the entire study period.Prior to analysis, the most dominant tree species was reclassified as one of 'ash, ' 'maple, ' 'cedar, ' 'other coniferous, ' or 'other deciduous' due to low numbers of observations for some deciduous and coniferous tree species as predominant among the sampling quadrats and to focus on effects of species with potential associations identified by Talbot et al. 9 .In all analyses, we mean-centred continuous variables: the proportion of canopy cover, soil moisture, and mean deer foraging intensity.
We performed all statistical analyses using the 'lme4' package in R version 4.1.3.Given the repeated transect measures performed within study neighbourhoods, we built generalized linear mixed models (GLMMs) with a random intercept for neighbourhoods to account for between-neighbourhood differences in the density and infection prevalence of ticks, as well as unmeasured differences in the configuration of neighbourhoods affecting relative distance between zones within neighbourhoods.In each model, we also forced month as an ordinal variable (5 to 10, May to October) to account for seasonal variation in peak activity of different tick life stages, and we used year as a factor with two levels (2020 and 2021) as an adjustment for interannual variation in tick abundance.To determine the most appropriate statistical distribution for our count data, we used the 'glmer' function to evaluate the association of our main predictor of interest, zone type, with tick densities in Poisson GLMMs.We checked these models for excess observations of zero ticks with the 'check_overdispersion' function from the 'performance' package.As overdispersion was detected, we proceeded with a negative binomial (NB) GLMM via the 'glmer.nb'function and checked the dispersion statistic again.Since the NB distribution adjusted for the overdispersion sufficiently, we deemed a zero-inflated process unnecessary.For models assessing associations with the presence of infected nymphs and adults, we used a Binomial GLMM with the 'glmer' function to estimate B. burgdorferi prevalence across our study sites.
For each outcome measure, we tested multivariable conceptual models conceived a priori.Negative binomial GLMMs exploring associations with the densities of nymphal and adult ticks included categorical variables (see Table 5 for ecological category definitions) for the dominant tree species, depth of leaf litter, understory density, and zone type, as well as continuous measurements for soil moisture, mean deer foraging intensity, and canopy cover.In Binomial GLMMs exploring the B. burgdorferi infection prevalence of nymphs and adults, we also included continuous predictors for the abundance and the B. burgdorferi infection prevalence of P. leucopus mice.We initially considered including predictors for P. maniculatus mice but omitted them due to their lower relative abundance between neighbourhoods and the near-total absence of B. burgdorferi infection among the sampled deer mice.A check of variance inflation factors among our model terms with the 'performance' package revealed collinearity between P. leucopus infection prevalence and zone type in both Binomial GLMMs.Thus, B. burgdorferi infection prevalence in white-footed mice was removed from the final models.For both count models, we calculated predicted tick densities based on the average marginal effects for neighbourhoods and

Figure 1 .
Figure 1.Location of residential (solid yellow fill), interface (light blue hatched-line fill), and woodland (bright green dot-hatched fill) sampling zones within each neighbourhood of western Ottawa: N1, N2, N3, and N4.The depicted zonal boundaries are approximations that include the locations of each drag sampling transect used during data collection, drawn to highlight their location and aid in visibility on this map.Basemap image is Copernicus Sentinel-2 data (2023), accessed and used in accordance with the regulations of the European Union.Context map base layer by Stamen Design, under CC BY 4.0, with data by OpenStreetMap, under the Open Database License (ODbL).

Figure 2 .
Figure 2. Average marginal effect on nymphal density (A), density of infected nymphs (B), tick (nymph and adult) density (C), and density of infected adults and nymphs (D) for interface and woodland zones compared to residential, with all other covariates at their mean value.Plots include the zonal effect within each neighbourhood (N1-N4) as well as independent of neighbourhood.The density of infected nymphs (B) and infected adults and nymphs combined (D) is determined by the product of the predicted density and the observed infection prevalence reported in Table1.

Table 1 .
58, Total Ixodes scapularis ticks collected of each life stage with mean nymphal and all-stage Borrelia burgdorferi-infection prevalence and 95% confidence intervals (95% CI) by neighbourhood and zone.

Table 3 .
Proportions (%) of sampling quadrats with the recorded observation during June 2020 and 2021 site visits within each neighbourhood and each zone type, with descriptive statistics for continuous measurements.

Table 4 .
Rate ratios (RR) for associations with I. scapularis nymphal and nymph and adult density and odds ratios (OR) for associations with B. burgdorferi infection prevalence of I. scapularis nymphs and ticks.All models are adjusted for sampling month, year, and repeated measurements in neighbourhoods.Estimates are reported with Wald 95% confidence intervals (95% CI) and P values.

Table 5 .
Ecological observations and measurements recorded at sampling locations, with description and rationale for consideration. )